Shear-Induced Stress Relaxation in a Two-Dimensional Wet Foam 
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We report on experimental measurements of the flow be- 
havior of a wet, two-dimensional foam under conditions of 
slow, steady shear. The initial response of the foam is elastic. 
Above the yield strain, the foam begins to flow. The flow con- 
sists of irregular intervals of elastic stretch followed by sudden 
reductions of the stress, i.e. stress drops. We report on the 
distribution of the stress drops as a function of the applied 
shear rate. We also comment on our results in the context of 
various two-dimensional models of foams. 
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Foams are ubiquitous in nature JJ and part of a larger 
class of materials which exhibit a type of behavior that 
has come to be called "jamming" Loosely speak- 

ing, a "jammed" material is one that is unable to flow, 
typically due to the packing of the constituent particles. 
For the case of foams, which are composed of gas bubbles 
separated by fluid walls, the jamming is a consequence 
of the topological constraints that develop as the bub- 
bles press against each other. The existence of a jammed 
state contributes to the complex flow behavior of foams. 
This flow behavior is one reason that foams are so inter- 
esting from both a fundamental and applied point of view 
[jlj. For sufficiently small strains or stresses, a foam acts 
as an elastic solid (the jammed state). However, when 
the strain (or the stress) exceeds a critical value, known 
as the yield strain, a foam begins to flow. Of particular 
interest is the flow behavior of foams for small rates of 
strain. This flow occurs through "avalanches" , or sudden 
nonlinear rearrangements of the bubbles in the foam that 
correspond to a decrease in the average energy, or stress, 
of the foam. In this paper, we will refer to these events 
as "stress drops" . This type of behavior is common to 
many jammed systems where flow occurs in an irregular 
or stick-slip type manner when the system is near the 
transition to jamming J^]. One of the outstanding ques- 
tions for flowing foams is the nature of the distribution 
of the sizes of the stress drops during this irregular flow. 

One reason that this remains an open issue is that 
predictions for the distribution of stress drop sizes are 



somewhat model dependent. Various models of two- 
dimensional foams have been used to study the connec- 
tion between topological rearrangements and flow prop- 
erties as a function of applied rate of strain |rj-[l^1. All 
of these simulations predict the same qualitative behav- 
ior for the macroscopic flow for sufficiently small rates of 
strain that was described above: elastic behavior below a 
yield strain and nonlinear bubble rearrangements above 
the yield strain. However, the quantitative predictions 
for the size distribution and frequency of these events are 
model dependent. The dependence is predominately due 
to different assumptions concerning the source of dissipa- 
tion in foams and treating foams with different degrees of 
"dryness". Distinguishing between these different mod- 
els, and their underlying assumptions, experimentally is 
an important step in our understanding of foam rheology. 

Two-dimensional foams may be characterized by the 
area fraction of gas, <j). For <fi < 0.84, foams "melt" into 
a froth of exclusively circular bubbles. Near this transi- 
tion, where bubbles are predominately circular, a foam 
is said to be "wet". In the limit <j) approaches 1, the 
bubbles become polygonal with infinitely thin walls, and 
a foam is said to be "dry". The simulations studied by 
Kawasaki, et al. [^-^] are based on the vertex models 
of foams. This model assumes infinitely thin walls and 
only treats the dynamics of the vertices. As such, it is 
applicable to dry foams and unlikely to apply to the wet- 
foam system studied in this work. The simulations of 
Weaire, et al. QJl^] focus on the behavior of foams un- 
der quasi-static, extensional flow. These are applicable 
to wet foams and, because they are quasistatic, contain 
no dissipation. These simulations measure the number 
of Tl events. A Tl event is a topological change that in- 
volves the switching of neighbors between bubbles. Ref- 
erence [[l2| reports events that include large numbers of 
Tl events, implying a distribution of stress drop with a 
power-law like behavior. In contrast, the bubble model 
of Durian, which is also applicable to a wet foam, pre- 
dicts power-law distributions with an exponential cutoff 
||[l0|. This model includes viscous dissipation between 
bubbles. Further work on this model demonstrated that 
the exponential cutoff increases as one approaches the 
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melting point |Q. However, the dependence on area 
fraction is weak, and true power-law behavior was never 
observed, even extremely close to the melting point [[L4|. 
Finally, simulations of a q-potts model of a sheared foam 
]l3| did not assume a particular form for the dissipation. 
This work suggests that the distribution of topological 
rearrangements is not power-law like; however, the distri- 
bution of energy drops may be consistent with power-law 



behavior |13|. 

In order to distinguish between these various the- 
oretical models, there has been limited experimental 
work on two-dimensional model foams ]l5| , p!6[ and three- 
dimensional foams ]l7|| . The work in Ref. |17| made in- 
direct measurements of the bubble rearrangements and 
concluded that only local rearrangements occurred dur- 
ing flow. The work in Ref. [|l5| used Langmuir monolayers 
to form two-dimensional foam. This work directly mea- 
sured the distribution of Tl events and concluded that 
the size distribution was not consistent with a power- 
law. In contrast, the work in Ref. @ used a single layer 
of bubbles trapped between a glass plate and a water 
surface. This work focused on the bubbles whose topo- 
logical class had changed. These measurements suggest 
that very large events can occur. 

The apparent discrepancies in these measurements 
may not be as severe as it first appears. Detailed simu- 
lations of various models show that the size of Tl events 
and the number of bubbles involved in a rearrangement 
are not necessarily correlated |l3|,[l4|]. Also, it is possible 
that the number of Tl events is not correlated with the 
size of a stress drop This is reasonable when one 

considers that not all Tl events will relieve an equivalent 
amount of stress. Therefore, to address the question of 
the scaling of the size distribution, it is critical to have 
direct measurements of the stress or energy drops. We 
have accomplished this for the case of another model, 
two-dimensional foam: a bubble raft fL8| , |l9f . 

A bubble raft consists of a single layer of bubbles 
placed on the surface of water. Bubble rafts have been 
used to model the flow behavior of amorphous solids 
[fL8lJl9| . The bubble rafts are an ideal system for the 
study of two-dimensional foams for a number of reasons. 
They allow for both direct measurements of the macro- 
scopic properties of the foam and the bubble dynamics. 
There are no confining glass plates, which can add com- 
plications as the shear rate is increased. Finally, one can 
control both the degree of order in the foam (by varying 
the distribution of bubble sizes) and the density of the 
foam with great precision. In this paper, we will report 
on the flow behavior of a disordered bubble raft with 
0^0.9. 

We generate flow of the bubble raft using a two- 
dimensional Couette viscometer that is described in de- 
tail in Ref. [^0| . The apparatus consists of two concentric 
cylinders oriented vertically. Water is placed between the 
cylinders, and the upper surface of the water is free. The 



outer cylinder consists of 12 individual pieces, so it can 
be expanded and compressed to adjust the density of the 
bubble raft. The working radius of the outer cylinder 
was 7.43 cm. The inner cylinder is two pieces: a solid 
cylinder placed in the water and a second piece that is 
in contact with the bubble raft and fits over the solid 
cylinder. The second piece is hung by a torsion wire and 
is used to measure the stress, a, on the inner rotor due 
to the bubble raft. It is important to note that this is a 
"two-dimensional stress" given by the force on the inner 
rotor in the tangential direction divided by the circum- 
ference of the rotor. Therefore, in terms of the torque, 
r on the rotor and the radius of the rotor r, the stress 
is given by a = r/(27rr 2 ). The radius of the inner cylin- 
der was 3.84 cm. A constant rate of strain is applied to 
the system by rotating the outer cylinder at a constant 
angular speed in the range 0.0005 rad/s to 0.01 rad/s 

The stress on the inner cylinder is determined from 
the angular displacement of the torsion wire support- 
ing the inner cylinder. The angular displacement was 
measured using magnetic flux. A coil was attached to 
the torsion wire and suspended within a high-frequency 
magnetic field. The induced voltage was used to deter- 
mine the angle of the coil. Typical values of the angle 
during shear ranged from 22° to 40°, corresponding to 
stresses in the range of 4 - 23 dyne/cm. The resolution 
in stress was set by our resolution in angular measure- 
ment and the torsion constant, n = 570 dyne cm/rad. 
The voltage signal was digitized using a 12-bit A to D 
converter in the computer. As a test of the noise level 
in the signal, the stress was monitored as a function of 
time without shear. In the absence of shear, the noise in 
the measured stress signal was at the level of the lowest 
bit in digitized signal, corresponding to changes in the 
stress of ±0.026 dyne/cm. Therefore, when measuring 
the changes in the stress, changes of ±0.026 dyne/cm 
were ignored, providing a lower limit on the size of the 
stress drops. 

The bubble raft is generated by flowing nitrogen 
through a solution of 44% by weight glycerine, 28% by 
weight water and 28% by weight Miracle Bubbles (Im- 
perial Toy Corp.). The bubble size is fixed by the pres- 
sure and needle diameter. For the experiments reported 
on here, three different size bubbles (2 mm, 3 mm, 5 
mm) were used, with approximately 500 bubbles of each 
size. The experiments were carried out over two hours. 
The coarsening of the foam without shear was monitored 
during this time, and no significant coarsening or stress 
drops due to coarsening were observed. However, after 
two-hours, the stability of the foam deteriorated rapidly, 
and shear resulted in severe rupture of bubbles. 

Figure 1 is an image of a section of the bubble raft 
between the two cylinders. The bubble raft was moni- 
tored for slippage at both the outer and inner cylinder. 
For all of the experiments reported here, there was no 
slip between the first row of bubbles and the correspond- 
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ing cylinder. Also, it was clear from the images that the 
stress drop events corresponded to rearrangements of the 
bubbles in the bulk of the foam. The average azimuthal 
velocity distribution was consistent with what one would 
expect for Couette flow. This is illustrated in Fig. 2 
for an angular rotation of 0.01 s _1 . The dashed line is 
the theoretical velocity profile for a Newtonian fluid, and 
the solid curve is the theoretical velocity profile for a 
shear-thinning fluid with a viscosity given by r] = ray 1-1 , 
where m is a constant, 7 is the shear rate, and n = 1/3 
pl[ . It should be noted that the bubble model predicts 
shear-thinning with an exponent of approximately 1/3 
|^2| . The fact that the velocity profile is consistent with 
the expected shear-thinning profile, and not the velocity 
profile of the water substrate, is strong evidence that the 
foam is acting independent of the water. Because the 
rate of strain varies as a function of radius, all reported 
values of for rates of strain are taken at the inner cylin- 
der. Finally, it should be noted that this is very different 
behavior from that observed for two-dimensional foams 
confined between glass plates. In this case, exponential 
decay of the velocity profile was observed ]23f| . 

Figure 3 shows a typical response of the bubble raft 
to shear. The stress on the inner cylinder is shown as 
a function of the applied strain for a rate of strain of 
3.1 x 10 -3 s _1 . The two key features to notice are the 
initial elastic region where the stress increases linearly 
with strain and the subsequent region in which flow oc- 
curs. The intermittent nature of the flow is obvious in 
this plot. The yield strain was always of order 1 for all 
rates of strain. 

Figure 4 shows the distribution of the stress drops for 
three different rates of strain. Here the size of the drop 
is normalized by the average stress per bubble. Also, as 
discussed above, any change in stress of ±0.026 dyne/cm 
was ignored and does not enter into the plotted distri- 
bution. The distribution of stress drops is not consistent 
with a power-law for all length scales; however, for small 
stress drops, the distribution is consistent with a power 
law with an exponent of -0.8. For comparison, the ex- 
ponent predicted for the bubble model is —0.70 ± 0.05 
lllof . Also, as predicted by the bubble model Q, the dis- 
tribution is essentially independent of the rate of strain 
for small stress drops, and the large stress drop cutoff is 
weakly dependent on the rate of strain. This is difficult 
to see from Fig. 4, but is made clear by considering the 
average stress drop size. For example, for a rate of strain 
of 0.031 s -1 , the average stress drop (normalized by the 
average stress per bubble) is 27.9; whereas, it is 47.0 for 
a rate of strain of 0.31 s -1 . 

We have conclusively demonstrated that the general 
features of the distribution of stress drops observed in 
the bubble raft are in agreement with the simulations re- 
ported in Ref. (To) and [Q. Specifically, the cutoff in 
the distribution for large stress drops is clear. Also, the 
viscosity of the bubble raft is consistent with a shear- 



thinning fluid with an exponent of 1/3. Finally, the dis- 
tribution of stress drops, particularly for the small stress- 
drops, is essentially independent of the rate of strain for 
the values studied here. This strongly suggest the exis- 
tence of a quasi-static limit. Such a limit is particularly 
important for our geometry where the rate of strain varies 
across the system. 

A number of important questions remain. First, what 
is the correspondence between the size of the stress drops 
and the number of bubbles that exchange neighbors? It 
is still possible that the larger stress drops are due to re- 
arrangements involving a wide distribution of the num- 
ber of bubbles, as reported in (l3| . Therefore, "system- 
wide" events are still a possibility. The current imaging 
system was not appropriate for careful tracking of bub- 
bles throughout the system. Future improvements in the 
imaging of this system will enable this question to be ad- 
dressed. Second, a detailed study of the flow behavior of 
bubble rafts as a function of both disorder (polydispertiy) 
and area fraction will provide further tests of the various 
two-dimensional models of foams. Of particular impor- 
tance is the behavior as a function of area fraction to 
distinguish further between the predictions of the quasi- 
static simulations QJl^] and the bubble model [ |To| , [l4[ . 
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FIG. 1. This is an image of one section of a typical bubble 
raft. Part of both the inner and outer cylinder are visible. 
The black scale-bar in the lower right corner is 3.6 mm. 





FIG. 3. Plot of the stress versus strain for a rate of strain 
of 3.1 x 10" 3 s _1 . 
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FIG. 4. Distribution of stress drops for three different rates 
of strain. Solid triangles are a rate of strain of 0.031 s _1 , solid 
squares are a rate of strain of 0.31 s _1 , and open circles are 
for a rate of strain of 0.48 s _1 . The solid line has a slope of 
-0.8 and is a guide to the eye. 
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FIG. 2. Plot of average azimuthal speed of the bubbles 
versus radial position (solid symbols). The solid line is the 
theoretical velocity profile for a shear-thinning fluid with n 
= 1/3 (see the text for details), and the dashed line is the 
theoretical profile for a Newtonian fluid. 
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